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Abstract 

We numerically study the existence of travelling breathers in Klein- 
Gordon chains, which consist of one-dimensional networks of nonlinear 
oscillators in an anharmonic on-site potential, linearly coupled to their 
nearest neighbors. Travelling breathers are spatially localized solutions 
having the property of being exactly translated by p sites along the chain 
after a fixed propagation time T (these solutions generalize the concept 
of solitary waves for which p = 1). In the case of even on-site potentials, 
the existence of small amplitude travelling breathers superposed on a small 
oscillatory tail has been proved recently (G. James and Y. Sire, to appear 
in Comm. Math. Phys., 2004), the tail being exponentially small with 
respect to the central oscillation size. In this paper we compute these so- 
lutions numerically and continue them into the large amplitude regime for 
different types of even potentials. We find that Klein-Gordon chains can 
support highly localized travelling breather solutions superposed on an os- 
cillatory tail. We provide examples where the tail can be made very small 
and is difficult to detect at the scale of central oscillations. In addition we 
numerically observe the existence of these solutions in the case of non even 
potentials. 

Keywords: travelling breathers; nonlinear lattices; numerical continuation, center 
manifold reduction. 
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1 Introduction 



We consider a chain of nonlinear oscillators described by the following system 
(discrete Klein-Gordon equation) 

cPx 

+ V'(x n ) = j(x n+1 + x n _i - 2x n ), n e Z. (1) 

System (0) describes an infinite chain of particles linearly coupled to their nearest 
neighbors, in the local potential V. We denote by x n the displacement of the nth 
particle from an equilibrium position, V(x) a smooth anharmonic on-site potential 
satisfying V"(0) = 0, ^"(0) = 1, and 7 > the coupling constant. This class 
of systems has been used in particular as a prototype for analyzing the effects 
of nonlinearity and spatial discreteness on localization of vibrational energy and 
energy transport (see [21] for a review). Equation (JJJ) yields in particular the 
Peyrard-Bishop model for DNA denaturation [35] ^] if one chooses V as a Morse 
potential. 

This paper deals with solutions of (JJJ) satisfying 

x n (r) = x n - p (r - T), (2) 

for a fixed integer p > 1 (p being the smallest possible) and Tel. The case 
when p = 1 in (J2J) corresponds to travelling waves with velocity 1/T. Solutions 
satisfying (J2j) for p ^ 1 consist of pulsating travelling waves, which are exactly 
translated by p sites after a fixed propagation time T and are allowed to oscillate 
as they propagate on the lattice. Solutions of type (J2J) having the additional 
property of spatial localization (x n (r) —>■ as n —>■ ±00) are known as exact 
travelling breathers (with velocity p/T) for p > 2 and solitary waves for p = 1. 

Approximate travelling breather solutions (i.e. spatially localized solutions 
satisfying (J2J) only approximately) have been extensively studied. These solutions 
have been numerically observed in various one-dimensional nonlinear lattices such 
as Klein-Gordon chains |TJ 113 El 13 EH , Fermi-Pasta-Ulam lattices 01 EH EH 
and the discrete nonlinear Schrodinger (DNLS) equation [TBI ESI ( see Ell for 
more references). These solutions have been studied also in more sophisticated 
nonlinear models describing energy transport in alpha-helix |2*^ and DNA 
[22] . One way of generating approximate travelling breathers consists of "kicking" 
static breathers consisting of spatially localized and time periodic oscillations 
(see the basic paper [3TJ for more details on these solutions). Static breathers 
are put into motion by perturbation in the direction of a pinning mode [HI E] . 
The possible existence of an energy barrier ("Peierls-Nabarro barrier") that the 
breather has to overcome in order to become mobile has drawn a lot of attention, 
see e.g. [TJI1 [21] and the review paper 00]. An analytical method recently 
developed by MacKay and Sepulchre [32] allows to approximate the nonuniform 
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breather motion along the chain (the breather center approximately moves under 
the action of an effective potential). An application to the Fermi-Pasta- Ulam 
lattice can be found in [29| . 

In the small amplitude regime, spatially localized solutions of ( f l ]) -(|2 ]) have been 
studied analytically by Iooss and Kirchgassner for p = 1 (case of travelling 
waves). Travelling wave solutions of (JTJ have the form x n (r) = Xq(t — nT) and 
are determined by the scalar advance-delay differential equation 



d 2 xo 



+ V'(x ) = 7 (x {t + T) - 2x + x (r - T)). 



(3) 



Extension to larger values of p has been performed by James and Sire |2H E2| 
(case of pulsating travelling waves). For fixed p > 1, problem (JH)-© reduces to 
the p-dimensional system of advance-delay differential equations 
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x 2 (t) - 2xi(r) + x p (t + T) 
x n +i(r) - 2x n (r) + x n _i(r) 



_ x x {t - T) - 2x p (r) + x p _i(r) _ 



(4) 



System Q is analyzed in the above references using a center manifold reduction 
in the infinite dimensional described e.g. in reference jlT]- It is rewritten 

as a reversible evolution problem in a suitable functional space, and considered 
for parameter values (T, 7) near critical curves where the imaginary part of the 
spectrum of the linearized operator consists of a pair of double eigenvalues and 
pairs of simple ones (the number of pairs is finite and equals the number of 
resonant phonons). Close to these curves, the pair of double eigenvalues splits 
into two pairs of hyperbolic eigenvalues with opposite nonzero real parts, which 
opens the possibility of finding solutions homoclinic to 0. Near these parameter 
values, the center manifold theorem reduces the problem locally to a reversible 
finite dimensional system of differential equations. The reduced system is put in a 
normal form which is integrable up to higher order terms. In some regions of the 
parameter space, the truncated normal form admits reversible orbits homoclinic to 
0, which bifurcate from the trivial state and correspond to approximate solutions 
of®. 

Exact small amplitude solutions of ((H)-(j2j) close to these approximate solutions 
(in a sense that we shall specify, allowing the existence of a small oscillatory tail) 
have been obtained for p = 1 and for p = 2 in the case of even potentials 
V [2E]- m both cases, the simplest homoclinic bifurcation yields a 6-dimensional 
reversible reduced system. As it is shown by Lombardi for different classes of 
reversible systems having a slow hyperbolic part and fast oscillatory modes (SDj 
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(see also jH] for related numerics), reversible solutions of the truncated normal 
form homoclinic to should not generically persist when higher order terms are 
taken into account in the normal form. The existence of corresponding travelling 
waves (p = 1) or travelling breathers (p = 2) decaying exactly to should be a 
codimension-1 phenomenon, the codimension depending on the number of pairs 
of simple purely imaginary eigenvalues in the considered parameter regime (there 
is one pair of purely imaginary eigenvalues, in addition to hyperbolic ones). How- 
ever, to confirm the nonexistence of reversible orbits homoclinic to (close to a 
small amplitude homoclinic orbit of the truncated normal form) for a given choice 
of V , 7, T, one has to check the nonvanishing of a certain Melnikov function be- 
ing extremely difficult to compute in practice [SO]. Due to this codimension-1 
character, in a given system (£Q) (with p = 1, or p = 2 and V even), 7 and V 
being fixed, exact travelling wave or travelling breather solutions decaying to 
at infinity might exist in the small amplitude regime, but for isolated values of 
the velocity p/T (see [20] for examples of nonlinear lattices supporting explicit 
travelling breather solutions with particular velocities). 

Instead of orbits homoclinic to 0, the full normal form admits orbits homo- 
clinic to small periodic ones jHO] (originating from the pair of purely imaginary 
eigenvalues). These solutions correspond to exact solitary wave {jp — 1) or trav- 
elling breather (p = 2) solutions of superposed on a small periodic oscillatory 
tail, which can be made exponentially small with respect to the central oscilla- 
tion size (the minimal tail size should be generically nonzero for a given value 
of (T, 7)). This phenomenon is in accordance with numerical observations in the 
case of travelling waves jH], and should explain numerical convergence problems 
arising in the computation of fully localized solitary waves or travelling breathers 
in various nonlinear lattices O |3] • 

For asymmetric potentials V and p > 2, the simplest homoclinic bifurcation 
yields a higher-dimensional (2p + 4-dimensional) reduced system, with supple- 
mentary pairs of simple imaginary eigenvalues of the linearized operator (the 
imaginary part of the spectrum consists of a pair of double eigenvalues and p 
pairs of simple ones). By analogy with results of Lombardi |30j . it has been 
conjectured j2Hl E2] that the above mentioned solutions of the truncated normal 
form homoclinic to do not generically persist when higher order terms are taken 
into account in the normal form. Persistence might be true if parameters (T, 7, 
coefficients in the Taylor expansion of V) are chosen on a discrete collection of 
codimension-p submanifolds of the parameter space. For general parameter val- 
ues, instead of orbits homoclinic to one can expect the existence of reversible 
orbits homoclinic to exponentially small p— dimensional tori, originating from the 
p additional pairs of simple purely imaginary eigenvalues. These solutions should 
constitute the principal part of exact travelling breather solutions of (Q) super- 
posed on a small quasi-periodic oscillatory tail. However, in order to obtain exact 
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solutions one has to prove the persistence of the corresponding homoclinic or- 
bits as higher order terms are taken into account in the normal form. This step 
is non-trivial and would require to generalize results of Lombardi jHOj available 
for a single pair of simple imaginary eigenvalues (another promising approach is 
developed in the recent work [26 J. 

One might think that this nonpersistence picture contrasts with the case of the 
integrable Ablowitz-Ladik lattice pQ, which supports an explicit family of exact 
travelling breather solutions in which the breather velocity can be varied contin- 
uously. However these solutions are not robust under various non-Hamiltonian 
reversible perturbations as shown in jH]. 

In this paper we numerically solve system (J1J in the case p — 2, which con- 
stitutes the simplest situation for travelling breather solutions (the case p = 1 
corresponding to travelling waves). This allows us to go beyond the analytical 
results obtained in 28\, which were restricted to the small amplitude regime. 

Our numerical observations confirm that system Q-(j2j) (with p = 2) admits 
travelling breather solutions superposed on a small oscillatory tail. These solu- 
tions are either weakly or strongly localized depending on the value of T (7 being 
fixed). These results are obtained both for hard and soft on-site potentials. In 
addition we provide examples where the tail's amplitude can be made extremely 
small near particular values of the breather velocity, which makes the tail difficult 
to detect at the scale of central oscillations. Note that travelling breather solu- 
tions have been previously computed in the DNLS equation (T7J |2] (in fact gauge 
symmetry reduces the problem to the computation of travelling waves), but the 
possible existence of a very small oscillatory tail has not been addressed. 

An intuitive explaination for the existence of a tail is that travelling breathers 
transfer energy to resonant phonons. The energy transfer can be balanced by the 
superposition of a resonant phonon tail, which allows the breathers to propagate 
in an exact way. This remark might suggest more generally that a way to en- 
hance breather mobility in nonlinear oscillator chains could be to excite resonant 
phonons with appropriate amplitudes. 

Our numerical computations are performed as follows. We consider system 
(HJ) for p = 2 and rescale time by setting t — We introduce the variable 
(ui(t), u 2 (t)) = {xi(t),x 2 (t+ -|)), which is a solution of the symmetrized system 



dt 2 



Ul 

u 2 



+ T 2 



v'( Ul ) 

V'(u 2 ) 



7T 2 



u 2 {t - \) - 2 Ul {t) + u 2 {t + \) 



ux{t + \)-2u 2 {t)+u x (t-\) 



V- 

21 J 



(5) 



and is related with the original variables by the relations 



Xti(t) = U\{^ — IL 2 —) if n is odd, 
%n(j) = w 2 (^ — if n is even. 



(6) 
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We solve system (j3J) by finite difference discretization and a Newton method (the 
Powell hybrid method 34 J. An alternative would be to use spectral methods, as 
done in [TH1 El EH Ej for the computation of travelling waves. 

In addition, we derive a simpler problem in the case of even potentials V. 
In this case, system © admits symmetric solutions satisfying ui = —u 2 , corre- 
sponding to solutions of (JTJ) having the form 



x*(t) = (-ir+V(^ - ^y^)- (7) 



This yields a new scalar advance-delay differential equation 



+ T 2 V'( Ul ) = -iT\ Ul (t + 1/2) + 2 Ul (t) + ui(t - 1/2)). 



The convergence of the Newton method heavily relies on having a good initial 
guess. For this purpose we use the explicit principal parts of solutions provided 
by center manifold theory in the small amplitude regime p8 j; and continue them 
numerically by varying T. 

An alternative numerical method would be to consider equation and search 
for fixed points of the time-T map F : (x n (0), x n (0)) i— > (x n+p (T), x n+p (T)), 
which correspond to solutions satisfying (j2j). This method has been used in jHj 
for p = 1 (case of travelling waves). The first step of the method consists in 
computing a static breather, which is then perturbed in the direction of a pinning 
mode in order to obtain a (slowly moving) approximate travelling breather. This 
solution serves in turn as an initial guess for the Newton method, as one computes 
fixed points of F. One limitation of the method is that pinning modes do not 
always exist (this is the case e.g. for certain hard quartic on-site potentials), but 
the method yields very precise results at least for travelling waves [S]. 

The paper is organized as follows. Section |21 describes the numerical scheme 
and suitable initial guesses for the Newton method. Section |3] is devoted to the 
scalar case (JHJ) • System (J3J) is numerically solved in section HJ 

2 Numerical method 

2.1 Finite difference scheme 

Equations (J3J) and (JHJ) have the following general form 

d 2 U 



dt 2 



+ W(U)= ^2AjU{t + J -), (9) 
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where U(t) £ R d (d = 1,2) and A-x,Ao,Ai are d x d matrix. The nonlinear 
function W : M d — > IR d is smooth and satisfies W(0) = 0. We consider periodic 
boundary conditions U(t + M) = U(t). For large M one can expect to find good 
approximations to spatially localized solutions which have "infinite period" (in 
practice we shall fix M = 20 most of the time). 

We shall approximate (jHJ) using a second order finite difference scheme. We 
choose the interval of discretization I — [— 4^ , 4^] and a discretization step h = 
where M, N £ N*. We fix an integer period M so that time delays in Q will be 
multiple of the step size, but this feature is not essential (see section . 

We set Ut = U(hi) with hi £ / and i = -MN + 1,...,MN. The finite 
difference scheme writes 



for i = —MN + 1, MN. We impose the following periodic boundary conditions 



Equations ()10 |) -(fTT j) lead to a nonlinear system of 2MN equations with 2MN 
unknowns t/j. We use a Powell hybrid method [33] to solve this system (we 
impose a stoping criteria at precision of 10~ 12 ). 

Numerical solutions of Q obtained by solving ([T0j) - fTT| are checked by an in- 
tegration of system ((H) with (x n (0), x n (0)) as initial condition (x n is deduced from 
equation (J0J) or (J2J), and one considers 2M lattice sites with periodic boundary 
conditions). The numerical integration is performed with a 4 th order Runge-Kutta 
scheme. We impose that x n is shifted by 2 sites after integration over [0, T] with 
a relative error tolerance of order 10~ 5 . Figure ^ gives a typical profile of the 
evolution of the relative error as a function of step h (note that the error decays 
as h 2 ). 

The crucial step is the choice of a good inital guess for the Powell method. 
This aspect is examined in next section. 

2.2 Initial guess 

The existence of small amplitude travelling breather solutions of (JIJ (having a 
small oscillatory tail) has been proved in [21] in the case of even potentials, using 
a center manifold technique. The method also provides leading order approximate 
expressions for these solutions (these approximate expressions are also available 
in the case of asymmetric potentials). This material is described below. 

We consider the linearization of system (jSJ) at {111,112) = (0,0) and search for 
solutions of the linearized problem in the form (u\, 112) = e tqt (iii, ii2) T ■ This leads 




(10) 



Ui + 2MN — Ui. 



(11) 
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Inh 



Figure 1: Relative error E = ||{x n (T) — x„_ 2 (0)}|| / ||{x n (T)}|| (|| || denotes the 
Euclidean norm) as a function of step h (in logarithmic scales), for a travelling 
breather solution at T = 7.7, 7 ~ 0.9, with the polynomial potential V(x) = 

1™2 _ 1™4 

to the following dispersion relation 



Equation (|T2"j) provides the spectrum (consisting in eigenvalues ±zg) of a certain 
linearized evolution operator obtained as system © is rewritten in the form of a 
reversible differential equation in a Banach space (see [2H] for more details). The 
spectrum on the imaginary axis corresponds to solutions q G E of (fT2"j) and is 
sketched in figure El 

Double roots gGMof (|T2")) satisfy in addition 



(solutions of (fT2" j) - (fT3"j) correspond to eigenvalues of the linearized evolution op- 
erator which are in general double, and at most triple [28 ). We note that for 



T 2 (l + 2 7 ) = (2k + 1) 2 tt 2 (k G N), q — (2k + 1)tt is a double root of JTJ). 



This situation occurs on the curves denoted as S 4 in figure 121 and corresponds 
to the crossing of two simple roots of (fT2*|) as these curves are crossed in the pa- 
rameter plane. The other double roots defined by system (|12 p -(|13 )) are found as 
parameters lie on the curve T : (T(q), r y(q)) (q G M + ) defined by 



(-q 2 + T 2 (1 + 2 7 )) 2 - 4( 7 T 2 ) 2 cos 2 (g/2) = 0. 



(12) 



2g(-g 2 + r 2 (l + 2 7 )) = (7T 2 ) 2 sin(g) 



(13) 



T 2 = q 2 -4gtan(g/4) 



(14) 



2q 



(15) 



7 



T 2 sin(?/2) 
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Figure 2: Sketch of the eigenvalues iq on the upper imaginary axis (right), in 
different regions of the parameter space (left). These eigenvalues correspond to 
real solutions of ()12j) (here q > 0). 



if q G [Akir, (2k + 1)271"] (for an integer k > 1), and 

T2 = o 2 + 1^/4)> < 16 > 

7 = ~T 2 sin(g/2)' ^ 

if q G [(2A; — l)27r,4A;7r] (k > 1). The range of q is determined by the condition 
T 2 > 0. One can observe that T lies in the parameter region 7T 2 > 4. 

We denote by the restriction of T to the interval q G [2/c7r, 2{k + l)7r] (these 
curves form "tongues" as shown in figure |2j). Our numerical computations will 
start near I\ where bifurcation of small amplitude travelling breathers take place. 
Travelling breathers bifurcate more generally near T 2 k+i j28j . and bifurcating 
solitary waves are found near Y 2 k [25J. 

Double roots on V correspond to the appearance (or vanishing) of a pair of 
real roots of (fT2"j) as the curve V is crossed. The envelope A of V is the bold 
line depicted in figure EJ Below it (and outside S 4 ), real solutions of |T2|) consist 
in two pairs of simple roots ±q±, ±q 2 . An additional pair of real (double) roots 
±qo appears as one reaches A from below (they correspond to pairs of hyperbolic 
eigenvalues of the linearized evolution operator colliding on the imaginary axis). 
In the sequel, we shall exclude points of A which are close to points where sqo + 
rqi + r'q 2 = for s,r,r f G Z and < \s\ + |r| + \r'\ < 4 (such values correspond 
to strong resonances), and denote this new set as A . 
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Now we fix (T ,7o) G T m n A . For (7, T) ~ (70, To), it has been proved in 
that small amplitude solutions of (J3J) have the form 

( 2$ ) = m ( ( ~i r ) +c{t) ( ) +jD(t) ( I ) + c - c -+*Mt),i,n 

(18) 

where u c = (A,B,C,D,A,B,C,D) T G C 8 and ^ : C 8 x R 2 -> M 2 is a smooth 
function, with *(w c ,7,T) = 0(||w c || 2 + ||« c ||(|7 - 7o| + |T — T |)) as (u c , 7,T) -> 
(0, 70, T ) (i.e. $ contain higher order terms). The complex coordinates (A, B, C, D) 
satisfy a reversible differential equation 

dA 

— = iq Q A + B + iAP(v 1 ,v 2 , v 3 , v,) + 0((\A\ + \B\ + |C| + |D|) 4 ), 
dB 

— = iq B + iBV(v 1 ,v 2 , v 3 , v 4 ) + AS(v u v 2 , v 3 , v A ) 
+0{{\A\ + \B\ + \C\ + \D\Y), 
■ iqiC + iCQ(v u v 2 , v 3 ,v 4 ) + 0((\A\ + \B\ + |C| + \D\) 4 ), (19) 

2g 2 TJ + iDT( Vl ,v 2 , v 3 , v 4 ) + 0{{\A\ + \B\ + |C| + |TJ|) 4 ), 



dC 

~dt 
dD 

~dt 



where 

vi = AA, v 2 = CC, v 3 = DD, v 4 = i(AB - AB) 

and V, S, Q, T are affine functions with smoothly parameters dependent real co- 
efficients, for (T, 7) in the neighborhood of A . We have 

V{vi,v 2 , v 3 , V4) = pi(7, T) + p 2 v x + p 3 v 2 + p 4 v 3 + p 5 v 4 , 

S(V1,V 2 , V 3 ,V 4 ) = Si(7,T) + S 2 V! + s 3 v 2 + s A v 3 + S 5 V4, (20) 
Q(ui, u 2 , u 3 , u 4 ) = ft (7, T) + g 2 t>i + g 3 t; 2 + g 4 w 3 + g> 4 , 

T(Mi, W 2 , «3, M4) = *l(7, r ) + *2«1 + ^2 + + hVA 

where pi,si,ft,£i vanish on Ao- The method for computing the coefficients of 
()20|) is given e.g. in [21], section 5.2 (see [21] f° r a description of the method in 
a general setting). Defining the coefficients in the expansion of V as 

V(x) = ^x 2 + |x 3 + ^x 4 + 0(|x| 5 ), x^O, (21) 

one has in particular for (T, 7) = (To, 70) G A 

go 4a 2 T 2 
(2 " = T ° 2( ~ 6/3 + " 27oT 2 cos(g ) - T 2 (l + 2 7o ) + A^' (22) 
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Moreover, 



as (T, 7) — > (T , 70) € A one has 



si = sn(7 - 7o) + Si 2 (T - T ) + h.o.t 



with 




e - cos(g /2) 



S12 = 4T 



(l + 27o)e + 27 cos(g /2) 
2e - ^ cos(g /2) 



(23) 



sn = - 



2e- 



^cos(g /2) 



We set e — — 1 in (J23~j) if qn F [(2fc - l)2?r, 4A;tt], jfe > 1 (i.e. (T ,7 ) G r 2fe _i) and 
e = 1 if g G [4A;7r, (2fc + l)2vr] ((T , 7 o) e T 2fc ). 



Equation (|I8j) expresses the fact that small solutions of (0) lie on a 8-dimensional 
center manifold -M 7i t> the dimension of which is equal to the number of marginal 
modes of equation linearized at (1*1,1*2) = for (7, T) = (70, T ). One can 
parametrize the center manifold in such a way that the reduced system (i.e. the 
8-dimensional differential equation satisfied by the coordinates u c of solutions on 
tM 7i t) is as simple as possible (this form is called "normal form"). The normal 
form of order 3 of the reduced system is given in (|19p. The truncated normal 
form (obtained by neglecting terms of orders 4 and higher in ([T9jl ) is integrable 



For (T, 7) w (T , 70) G T m f] A and if (T, 7) lies below A (this corresponds to 
having Si > in (j2Uj) ). the linearized operator at the right side of (fT9j) admits two 
pairs of symmetric hyperbolic eigenvalues close to ±2g > an d two pairs of purely 
imaginary eigenvalues close to ±iqi, ±2g 2 - If ■S2(7o,T ) < 0, the truncated normal 
form admits reversible homoclinic solutions to which bifurcate from the trivial 
state as (T, 7) — > (To, 70). They correspond near T 2 k+\ to approximate travelling 
breather solutions of (fl)). given by system (J5j) and having the form 



r (t) = (^i) 1 /2( cosh(ts i/2 )) -i ! ^ = pit + 2 ^2 s i/2 tanh ( ts V 2 ). (25) 



The ansatz (j2H) provides a good inital guess for computing small amplitude 
travelling breather solutions with the Powell method. The solutions initially 
computed are then followed by continuation as a function of T towards higher 
amplitudes. 



[23, m 





(24) 



and 
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The normal form coefficients pi,P2 in ()25p can be computed as explained 
in [21]. However, the term ip(t) in (J2"2j) is negligible as (T, 7) (To, 70) since 
Pi, si ~ and t G [— y> y]- Consequently one can fix ?/> = in (}24|) to initiate 
the Powell method. 

In addition, it is not necessary to take oscillatory modes C(t),D(t) into ac- 
count in the initial ansatz ()24|) because those are conjectured to be exponentially 
small with respect to (7 — 70 1 + |T — T | as (7, T) — > (70, To), for some exact 
solutions of (jHJ) (this has been proved for even potentials |28j). Note in addition 
that a tail appears naturally as one iterates the Powell method. 

In the following section, we present numerical computations in a particular 
symmetric case (for even potentials V) in which travelling breather solutions are 
given by a scalar advance-delay differential equation. 

3 Case of even potentials V 

We consider the case when the potential V in (JTJ is even. In this case, there exist 
solutions of (JTJ) satisfying 

T 

x n +i{r) = -x n (r - -), (26) 

given by the scalar advance-delay differential equation (|SJ). Section f3.1l reviews 
analytical results [2E] concerning small amplitude travelling breather solutions 
of the form (|26|) . These solutions are numerically computed in section 13.21 and 
continued into the large amplitude regime. In section I3~3l we numerically identify 
these solutions as modulations of certain time-periodic oscillations (time-periodic 
pulsating travelling waves). 

3.1 Exact small amplitude travelling breathers 

In this section we sum up analytical results obtained in ,28j, which follow from 
the center manifold reduction theorem and the analysis of the reduced equation 

The homoclinic orbits of the truncated normal form described in section 12.21 
bifurcate near parts of the "tongues" T m where S2 < 0. For even potentials (a = 
in ([21))). homoclinic bifurcations occur near left branches of T m if V is hard {(3 > 
in (|2*T|)). and near right branches if V is soft (/3 < 0). The situation is depicted 
in figure H3 (see [28 for more details). 

Now assume £2(70, T ) < defined by equation (|2^|) for a fixed (To, 70) G 
Ao P| T 2 k+i and consider (7, T) w (70, T ), (T, 7) lying below A in the parameter 
plane. The full reduced equation (fT^j) admits an invariant subspace defined by 
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Figure 3: Regions in the parameter space where small amplitude homoclinic orbits 
to bifurcate for the truncated normal form (case of an even potential). 

D = due to the evenness of V (see |2E])- Considering (j!9|) on this invariant 
subspace is equivalent to searching for solutions of (HI)-© satisfying (m = 
2k + 1 and u 2 = -Ui in (ffBjl ). 

If (T ,7o) lies outside some subset of Aon^it+i having zero Lebesgue mea- 
sure (corresponding to resonant cases), the full reduced equation (fT9|) admits 
small amplitude reversible solutions (with D = 0) homoclinic to periodic orbits. 
These solutions correspond to exact travelling breather solutions of system (JTJ 
superposed at infinity on an oscillatory (periodic) tail. Their principal part is 
given by 

T 11 — 1 

x n (r) = A + C] (I - — ) + c.c. + h.o.t, (27) 

where (A, C) is a solution of (|T9|) homoclinic to a periodic orbit (see [2Hj for an 
approximation of (A,C) at leading order). For a fixed value of (j,T) (and up 
to a time shift), these solutions occur in a one-parameter family parametrized by 
the amplitude of oscillations at infinity. The lower bound of these amplitudes is 
0(e~ c ^ 1/2 ), where fj, = \T - T \ + \j - j \, c> 0. 

The lower bound of the amplitudes should be generically nonzero, but may 
vanish on a discrete collection of curves in the parameter plane (T, 7) . As a 
consequence, in a given system (JTJ) (with fixed coupling constant 7 and symmetric 
on-site potential V), exact travelling breather solutions decaying to at infinity 
(and satisfying (J2J for p = 2) may exist in the small amplitude regime, for isolated 
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Figure 4: Large amplitude solution in the case of the trigonometric potential 
V{x) = 1 - cos(x) with T = 8.1,7 = 70 « 0.9, T ta 6.63 ((T ,7 ) G The left 
figure shows the displacements of mass n = 1 as a function of i (note that t is 
a rescaled time). The right one shows displacements for all lattice sites at time 
r = 0. 



values of the breather velocity 2/T. 

In the following section, we fix (To, 70) G r\ and consider \i = |T — T | as 
a small parameter, 7 = 70 being fixed. We numerically solve equation (JSJ) and 
follow the above mentioned travelling breather solutions into the large amplitude 
regime. 



3.2 Numerical computation of travelling breathers 

Let us consider the trigonometric potential V(x) = 1 — cos(a;) {(3 = — | in (j2U)) 
and fix 7 « 0.9, T = 8.1. These parameter values lie near the existence domain of 
small amplitude travelling breather solutions (see section l3~T|) . However, param- 
eters are chosen sufficiently far from the curve F± on which travelling breathers 
bifurcate (the point (T ,7 ) (6.63,0.9) lies on Ti), and correspond thereby to 
highly localized solutions. Figure |U shows a solution Ui(t) of (JSJ) (left), and the 
displacements of lattice sites at time r = deduced from equation © (right). 

This result shows that the Klein-Gordon lattice with a trigonometric potential 
supports highly localized travelling breather solutions. Although the solution 
in figure 0] seems perfectly localized, it admits in fact a small oscillatory tail 
not visible at the figure scale (this phenomenon is known for small amplitude 
solutions, see section Q). The tail is clearly visible after magnification (figure 

In the sequel we study in more detail the properties of these solutions as T 
is varied, and depending whether one considers a hard or soft on-site potential. 
For this purpose we shall fix V(x) = \x 2 + |x 4 . Figure |H1 presents numerically 
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Figure 5: Magnification of the solution tail in figure 0] (case of a trigonometric 
potential V,T = 8.1, 7 « 0.9). 

computed solutions of (JHJ) for different values of the propagation time T. Left 
column corresponds to the hard potential V(x) = \x 2 + ^x 4 (/3 = 1/6), the right 
one to the soft potential V(x) = \x 2 — \x A {(3 = —1). We start close to parameter 
values (T , 70) € I\ fl A and fix 7 = 70. 

One can check that as /z = \T — Tq\ increases, the central peak of the solution 
gets more localized and increases in amplitude. In addition, a periodic tail be- 
comes more visible as /1 increases. One can observe a frequency difference between 
the central part of solutions and their tail, as it is the case in the small amplitude 
regime (in equation (|27j). the central part A(t) and the tail C(t) oscillate respec- 
tively with frequencies close to the linear frequencies go and qi). In particular the 
tail has a lower frequency in the case of hard potentials, and a larger one for soft 
potentials. In figure [7| we plot the Fourier spectrum of u\(t) for two profiles of 
figure El (hard and soft potential case). The tail frequency and the central part 
internal frequency are found relatively close to the linear frequencies qi and go- 

In figure |H1 we compare the tail of a localized solution of figure |H1 with a 
numerically computed periodic solution (a nonlinear normal mode) originating 
from the eigenvalues ±iq% (the existence of such solutions is guaranteed by the 
Lyapunov-Devaney theorem). The tail of the localized solution is very close to 
the periodic solution, which confirms that the large amplitude localized solution 
can be seen as "homoclinic to periodic", as in the small amplitude regime. 

In figure H3 we plot as a function of n a travelling breather solution of 
corresponding (via equation (j7J)) to a profile of figure El In the case of a hard 
potential one observes that most neighboring oscillators are out of phase near 
the breather center. The case of a soft potential is considered in figure which 
shows that neighboring oscillators are in phase near the center. 

Now let us examine the amplitude of travelling breathers versus parameter \x 
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Figure 6: Localized numerical solutions of (jEJ) for various propagation times T 
and polynomial potentials V(x) = \x 2 + |x 4 . The left figures correspond to a 
hard potential (fi = 1/6, 7 = 70 w 0.83, T w 5.59 and T = 5.53, T = 5.5, 
T = 5.35 from top to bottom) and the right ones to a soft potential (j3 = — 1, 
7 = 70 w 0.9, T w 6.63 and T = 6.8, T = 7.4, T = 8.1 from top to bottom). The 
value of n — \T — T \ increases from top to bottom. 
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Figure 7: Fourier spectra of u\{t) (profiles of figure |B1 are extended with periodic 
boundary conditions). The right figure corresponds to a solution at T = 7.4,7 = 
7o ~ 0.9, T ~ 6.63 for V(x) = \x 2 — \x 4 (see figure El right). The left one 
corresponds to a solution at T = 5.5, 7 = 70 ~ 0.83, T ~ 5.59 for V{x) = \x 2 + 
^x 4 (see figure IHl left). The peaks are relatively close to qo/2-n rs 1.11, qi/2n ^ 
2.30 in the soft potential case and go/2vr m 1.7, gi/27r w 0.9 in the hard case. 




t 



Figure 8: Comparison of a nonlinear normal mode issued from the eigenvalues 
±iqi (continuous line) with the tail of a localized solution of figure El (dotted line), 
for T = 5.35 and V(x) = \x 2 + ±x A . We have 7 = 70 « 0.83, T w 5.59 and 
C/|,.'u) i- r,. 
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Figure 9: Numerical solution of CQ)-© for T = 5.5 and p = 2. We consider 
V(x) = \x 2 + ^x 4 , 7 = 70 « 0.83, T w 5.59 ((T ,7 ) G Tx). The solution 
is plotted as a function of n, for r = 0, and corresponds via equation (JJJ) to a 
solution U\{t) in figure (left column, middle). 




Figure 10: Numerical solution of (JH)-(j2J for T = 7.4 and p = 2. We consider 
V(x) = \x 2 -\x A } 7 = 70 w 0.9, T w 6.63 ((T ,7 ) G Ti). The solution is plotted 
as a function of n, for r = 0, and corresponds via equation © to a solution 
in figure |H1 (right column, middle). 
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Figure 11: Maximal amplitude HuiHoo °f a numerical solution of (jHJ) (continuous 
line) as a function of ^ in logarithmic scales. The dotted line represents a linear 
regression. One has V(x) = y + |j, 7 = 70 ~ 0.83, T ps 5.59 for the left figure 

2 4 

and V(x) — \ — j) 7 = 7o ~ 0.9, To ~ 6.63 for the right one. In both cases one 
has fj, — \T — T 1 and (T , 70) G IY 

in the large amplitude regime. In figure ^2 we plot the amplitude of solutions for 
polynomial potentials V(x) = y + jx 4 . One finds that the amplitude behaves 
approximately as // for relatively large values of \i (corresponding to highly local- 
ized solutions), with 5 ~ 0.44 for the hard potential (with f3 = 1/6) and 5 ~ 0.46 
for the soft one (with (5 = — 1). In this latter case the amplitude graph becomes 
less regular as one enters the next "tongue" T2 (figure |2J) , which might indicate 
bifurcations occuring near the numerically computed solutions (this part is not 
visible in figure ITTj) . 

A /i 1//2 dependency of the amplitude (at leading order) is expected in the small 
amplitude regime fi m |28j . Surprisingly a rather similar behaviour is observed 
in the large amplitude regime. 

We have found in addition multibump solutions. Figure ^] shows e.g. a 
2— bumps solution of equation (JHJ) (note the reverse symmetry between the two 
bumps) . 

3.3 Modulated periodic solutions 

As previously we consider (T , 70) G Ti n A and choose \i — |T — T | as a small 
parameter, 7 = 70 being fixed. If one neglects higher order terms in the normal 
form (|19|). the reversible bifurcation near Y2k+\ creates (in addition to homoclinic 
solutions) a three-parameter family of periodic solutions with C = D = and 

A, xe {t) = re^ t+d \ B^t) = i\re^ t+d \ (28) 



19 




Figure 12: 2— bumps solution of equation ©. We consider V(x) = \x 2 — jX 4 , 
T « 6.63, 7 = 7o « 0.9 ((T , 7 o) e i\) and T = 9. 

provided one approaches A from below in figure |21(si > 0), s 2 < and expression 
(|28|) is considered for A ~ (A is a real parameter). The frequency Q and modulus 
r depend on (/x, A) as follows 

n - q = A + Pl - (pa + 2Ap 5 ) A + ^ , (29) 

s 2 + 2As 5 

r 2 = - A2 + ^ . (30) 
s 2 + 2As 5 

Note that coefficients pi, Si depend on /x and the first integral v± of the truncated 
normal form takes the value V4 = i(A^\ t gB^\ t g — A^\ t gB^g) = 2r 2 A. The 
periodic solutions (|2*Hj) are 0{J~ji + |A|) as (/x, A) — > (si vanishes on A) and 
their frequency is close to go (the colliding pair of double eigenvalues is ±iqo) 
since pi vanishes on A. 

These solutions should correspond to solutions of (JTJ having the form 

T 11 — 1 

x n (r) « 2(-l)» r cos - ^-) + 0] (31) 

at leading order, and consisting in time-periodic pulsating travelling waves. One 
has equivalently 

(^))»2 r cos[« + 9 ](- 1 ). (32) 

To obtain exact solutions of (JTJ instead of leading order approximations, one 
should prove the persistence of the periodic solutions (|28j) as higher order terms 
are taken into account in the normal form (|19j). We shall not treat this question 
in the present paper and leave it for future works. 
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For computing the corresponding solutions of © with the same approach as 
in section 13.21 one has to consider the period M of solutions as a real parameter. 
The ansatz (J32|) is used with Q = kjr (for an integer k > 1) to initiate the 
Powell method. Equations (}29|) and (J30j) determine A,r k locally as functions 
of the two parameters /J, ~ and M « fc— (hence Q ~ g )- If we fix i n addition 
M — k 2 ^ = 0(\fj,\), then equation (|2^j) yields A = and we obtain r 2 = 

—S\/s2 + 0(fi 2 ) using equation (J30j) . Consequently we fix r 2 = — S1/S2 in the 
ansatz (J3~2~j) to initiate the Powell method near the "tongue" Ti. 

To adapt the numerical scheme f|T0^) - f|l If) to non-integer values of M, it suffices 
to rescale equation (JHJ) into 

+ M 2 TV(?/i) = -M 2 7 T 2 ( 2/1 (x + JL) + 2y 1 {x) + y x (x- JL)), (33) 

where y\{x) = Ui(t), x = and (|3*3*|) is subject to periodic boundary conditions 
yi(x+ 1) = yi(x). One discretizes (J33|) with a scheme similar to (fT0|) - fTT| . except 
advance and delay terms are now obtained by linear interpolation. 

We have numerically observed that localized solutions bifurcate from the above 
family of periodic solutions as M is fixed and T is varied (figure IT^j) . This property 
could be used to compute localized solutions by continuation from periodic solu- 
tions. This method has been previously employed in ^7] for computing travelling 
breathers in the DNLS system. In this case the problem reduces to the com- 
putation of solitary waves, which bifurcate from families of (explicitly known) 
periodic travelling waves. This approach seems less efficient in our case since 
periodic solutions are not explicitly known and have to be computed numerically. 

In the next section, we present numerical computations in the case of a non- 
even potential, corresponding to system (JBJ). 



4 Case of non even potentials V 
4.1 Approximate solutions 

The homoclinic orbits of the truncated normal form described in section 12.21 bi- 
furcate near parts of the "tongues" T m where s 2 < (see equation (J22J)- in the 
case on non even potentials, these regions vary with the parameter 

k = f3/a 2 } (34) 

depending on the quartic and squared cubic coefficients in the potential V (see 
(JHJ). The situation is depicted in figure ITU for the curve Ti (see |2Hj for the 
detailed study of the sign of s 2 ). 
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Figure 13: Comparison of a periodic solution (continuous line) with a weakly 
localized one (dotted line) for the soft polynomial potential V(x) = \x 2 — \x A 
and different values of T : T = 6.75 (top), T = 6.8 (middle), T = 7.2 (bottom). 
We have 7 = 70 « 0.9 and T w 6.63, (T , 70) G I\. 
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Figure 14: Part of I\ where S2 < (bold line). The dashed area is a small 
neighborhood of this curve lying below A. It corresponds to regions where small 
amplitude homoclinic orbits to bifurcate for the normal form (|19|) without higher 
order terms. The parameter k defined by (}34"J) varies as follows : (i) k < 10/9, 
(ii) 10/9 < k < k c with k c w 1.15, (iii) k c < k < 4/3, (iv) k > 4/3. 
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Now let us assume £2(70, Tq) < for a fixed (T , 70) G A P| r 2 £;+i and consider 
(7, T) « (70, T ), (T, 7) lying below A in the parameter plane. This is the case 
e.g. in the dashed area near Ti depicted in figure ITU 

Since V is non even, D = is not in general an invariant subspace of the 
full normal form (fT^j) and the oscillatory mode D{t) cannot be eliminated as 
in section 13.11 In fact, one expects (this is still a conjecture) that bifurcating 
localized solutions have in general nonvanishing (up to an exponentially small 
size) oscillatory tails D(t), C(t) oscillating at (fast) frequencies close to q 2 , q\- 
These oscillations are combined with fast oscillations at a frequency close to go 
(oscillatory part of A(t)) and the slow hyperbolic dynamics of the envelope of 
A(t) (time scales are of order 1/ ' J~ji, fi = |T — T | + (7 — 7 |). 

The normal form (|T9*j) truncated at order 4 admits small amplitude solutions 
homoclinic to 2-tori (originating from the pairs of eigenvalues Hqi, rkiq^). These 
approximate solutions of (fT9j) have an explicit form (see j2EJ, section 5.3). It 
has been conjectured in j2H] that they should correspond to the principal part of 
travelling breather solutions of system (0), superposed at infinity on an oscillatory 
(quasiperiodic) tail, and given at leading order by the expression 

T 71 — 1 

x n (r) « [ (-1)" A + (-1)" C + £>](- - — ) + c.c. (35) 

((A,B,C,D) denotes one of the above mentioned approximate solutions of ()19|) 
homoclinic to a small quasiperiodic orbit). The persistence of these homoclinic 
solutions for the full normal form is still an open mathematical problem, which 
would require to generalize results of Lombardi 30J available when one of the 
oscillatory modes C(t) or D(t) is absent. 

4.2 Numerical computation of travelling breathers 

In this section we consider the Morse potential V(x) = |(1 — e~ x ) 2 and perform 
the same type of computations as in section l3~2*l except we now work with system 
(jSJ) instead of the scalar equation (jHJ). We consider a point (To, 70) ~ (6.63,0.9) 
on Ti, fix 7 = 70 and let 11 = T — T vary. 

To initiate the Powell method we use the approximate travelling breather 
solutions described in section 12721 For this purpose one has to choose (T, 7) near 
parts of Fx where s 2 < 0, as described in section |4~H figure ITU ( dashed area). The 
Morse potential has the expansion (|2*Tj) with a = —3/2, (3 = 7/6, and k = 14/27. 

Figure [T51 presents numerically computed solutions of (0) for different values 
of T. We only plot U\(t) because U2{t) behaves similarly. As ji — T — T increases, 
the central hump of the solution gets more localized and increases in amplitude. 

Note that U\{t) and —U2{t) are different, contrarily to the case of even po- 
tentials. This asymmetry can be observed in figure It corresponds to the 
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Figure 15: Localized numerical solutions of (jSJ) for various propagation times 
T and the Morse potential V(x) = \{\ — e~ x ) 2 . We plot u\{t) {u2{t) behaves 
similarly). One has T = 6.7, T = 7.15 and T = 7.45 from top to bottom, 
(T , 70) ~ (6.63, 0.9) e Y\, 7 = 70. The value of fi = T — T increases from top to 
bottom. 
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Figure 16: Localized numerical solution (ui(t) , u 2 (t)) of (JSJ) for T = 7.45 and the 
(asymmetric) Morse potential V(x) = \{l—e~ x ) 2 . One has (T , 70) ~ (6.63,0.9) e 
Ti, 7 = 7o- Note the difference between U\(t) and —U2{t). 

fact that the symmetry x n+ i{r) = —x n {r — -|) is broken in the case of non even 
potentials. 

In figure El we plot as a function of n the corresponding mass displacements 
x n (r), at times r = and r = T (we fix T = 7.15 in the left figure and T = 7.45 
in the right one). As prescribed by condition (j2j), one recovers the same profile at 
both times, up to translation by 2 sites. One can note that neighboring oscillators 
are in phase near the breather center. 

At T = 7.15 (figures HS1 middle, and figure IT71 left) no tail is visible at least at 
the figure scale. On the contrary a tail appears for T = 7.45 (figures E3 bottom, 
and figure IT71 right). 

Figure provides a magnification of the tail of the solution profile. The 
oscillations observed near the boundary for T = 6.9 are not part of the solution 
tail and correspond in fact to oscillations of the central hump at a frequency close 
to go (near this value of T the solution is weakly localized, see figure ITBI top). The 
magnification shows small oscillations for T = 7.15 with a frequency close to the 
frequency observed for T = 6.9, therefore we attribute these oscillations to the 
damped central part of the solution, i.e. to the component A(t) approximated by 
(|24|). Consequently we do not find for T = 7.15 a visible track of an oscillatory 
tail corresponding to components C(t), D(t). On the contrary, in figure fTHl a 
nondecaying tail with a frequency close to q 2 is clearly visible at T = 7.45 and 
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Figure 17: Numerical solution of ©-© for T = 7.15 (left), T = 7.45 (right) and 
p = 2. We consider V(x) = |(1 - e^) 2 , 7 = 7o « 0.9, T « 6.63 ((T ,7o) e r^. 
The solution is plotted as a function of n, for r = (continuous line) and r = T 
(dotted line). It corresponds via system (j3J) to profiles tti(t) in figure fTBI f middle 
and bottom). 

larger values of T (i.e. for larger values of //). The nondecaying tail is also visible 
for some smaller values of T but the amplitude becomes very small (the tail size 
is decreased approximately by a factor of 100 for T = 7.35, see figure [THJ). 

FigurelT9lshows the Fourier spectrum of u\(t) for T = 7.45 and one can observe 
a sharp peak at a frequency close to ^2/2^ « 1.6. The damped oscillations of the 
central hump (approximated by (J2H0 at a frequency close to go appear as broad 
peaks and involve harmonics of all orders. In particular, figure ITBI ( bottom) clearly 
shows that the oscillatory part of the central hump has a nonzero time average. 
This situation contrasts with the even potential case (figure |7J), in which these 
oscillations only involve odd order harmonics. 

Another difference with respect to the even potential case is that solutions 
tails correspond here to the component D of (|33j) instead of C (the tail oscillates 
at a frequency close to qi instead of qi). Consequently the oscillatory tail consists 
in a time-periodic travelling wave, instead of a time-periodic pulsating travelling 
wave as in ()27|) . Note that we conjecture in fact the existence of a very small 
oscillatory component C(t) (not visible in figure IT6*j) . which might be hidden in the 
broad peak of figure EH close to 2q (|p ~ 2.22 and ^- ~ 2.30 are rather close). 
This is by analogy to the small amplitude regime described by the normal form 
(|19|). Indeed we conjecture that bifurcating homoclinic orbits of the full normal 
form have generically nonvanishing components C,D, except if T, 7 satisfy two 
(one for each component) compatibility conditions (|2Hj, section 5.3). 

We end by a study of travelling breathers amplitude versus /i = T — T in the 
large amplitude regime ffigurel2T)j). One finds that the amplitude behaves approx- 
imately as fi s for relatively large values of [i (corresponding to highly localized 
solutions), with S ~ 0.56. A jj 1 ^ 2 dependency of the amplitude (at leading order) 
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Figure 18: Magnification of the tail of localized solutions (component Ui(t)) for 
the Morse potential with 7 = 70 ~ 0.9, T ~ 6.63 ((T ,7 ) G Tx). We consider 
several values of T : T = 6.9 (crosses x), T = 7.15 (dashed line), T = 7.35 
(continuous line), T = 7.45 (crosses +). The minimal tail size is found around 
T = 7.15. Cases T = 7.15 and T = 7.45 can be compared with figure [TBI 




Figure 19: Fourier spectrum of Ui(t) (the spectrum of U2(t) is similar) for T = 
7.45, 7 = 70, (T , 70) ~ (6.63, 0.9) G T±. The corresponding profile of figure EH] 
has been extended with periodic boundary conditions. Note the sharp frequency 
peak close to f- w 1.6. One has f&. w 1.11, f- w 2.30. 
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Figure 20: Maximal amplitude ||wi||oo of a numerical solution of (J5J) (continuous 
line) as a function of /z, in logarithmic scales. The dotted line represents a linear 
regression. One has V(x) = \{l - e~ x f , 7 = 70 « 0.9, T « 6.63 ((T ,7 ) e T^, 
fi = T — T . A similar behaviour is observed for U2{t). 

is expected in the small amplitude regime /1 ~ |2H] • Surprisingly a rather similar 
behaviour is observed in the large amplitude regime. In addition the scaling in 
amplitude is rather similar to the case of equation (JHJ) with a quartic interaction 
potential (section l3~2|) . 

References 

[1] M.J. Ablowitz and J.F. Ladik. Nonlinear differential- difference equations and 
Fourier analysis, J. Math. Phys. 17, 1011-1018 (1976). 

[2] M.J. Ablowitz, Z. Musslimani and G. Biondini. Methods for discrete solitons 
in nonlinear lattices, Phys. Rev. E 65 (2002). 

[3] M.J. Ablowitz and Z. Musslimani. Discrete spatial solitons in a diffraction- 
managed nonlinear waveguide array : a unified approach, Physica D 184, 
276-303 (2003). 

[4] A. A. Aigner, A.R. Champneys and V.M. Rothos. A new barrier to the exis- 
tence of moving kinks in Frenkel-Kontorova lattices, Physica D 186, 148-170 
(2003). 



29 



[5] S. Aubry, T. Cretegny. Mobility and reactivity of discrete breathers, Physica 
D 119, 34-46 (1998). 

[6] A. Berger, R.S. MacKay, V.M. Rothos. A criterion for non-persistence of 
travelling breathers for perturbations of the Ablowitz-Ladik lattice. Discrete 
Cont. Dyn. Sys. B 4, n. 4, 911-920 (2004). 

[7] S.R. Bickham, S.A. Kiselev and A.J. Sievers. Stationary and moving intrin- 
sic localized modes in one-dimensional monoatomic lattices with cubic and 
quartic anharmonicity, Phys. Rev. B 47, 21 (1993). 

[8] A.R. Champ neys. Codimension-one persistence beyond all orders of homo- 
clinic orbits to singular saddle centres in reversible systems, Nonlineahty 14, 
87-112 (2000). 

[9] D. Chen, S. Aubry, CP. Tsironis. Breather mobility in discrete 4 lattices. 
Phys. Rev. Lett. 77, 4776 (1996). 

[10] J Cuevas, J.F.R. Archilla, Yu.B. Gaididei, F.R. Romero. Moving breathers in 
a DNA model with competing short- and long- range dispersive interactions, 
Physica D 163, 106-126 (2002). 

[11] J Cuevas, F. Palmero, J.F.R. Archilla, F.R. Romero. Moving discrete 
breathers in a Klein-Gordon chain with an impurity, J. Phys. A : Math. 
Gen. 35, 10519-10530 (2002). 

[12] T. Dauxois, M. Peyrard and A.R. Bishop. Dynamics and thermodynamics of 
a nonlinear model for DNA denaturation, Phys Rev E 47, n. 1, 684 (1993). 

[13] T. Dauxois, M. Peyrard and C.R. Willis. Discreteness effects on the formation 
and propagation of breathers in nonlinear Klein-Gordon equations, Phys. 
Rev. E48, 4768 (1993). 

[14] D.B. Duncan, J.C. Eilbeck, H. Feddersen and J.A.D. Wattis. Solitons in 
lattices, Physica D 68, 1-11 (1993). 

[15] J.C. Eilbeck. Numerical simulations of the dynamics of polypeptide chains 
and proteins. In Chikao Kawabata and A.R. Bishop, editors, Computer Anal- 
ysis for Life Science - Progress and Challenges in Biological and Synthetic 
Polymer Research, pages 12-21, Tokyo, 1986, Ohmsha. 

[16] J.C. Eilbeck and R. Flesch. Calculation of families of solitary waves on dis- 
crete lattices, Physics Letters A 149, 200-202 (1990). 



30 



[17] H. Feddersen. Solitary wave solutions to the discrete nonlinear Shrodinger 
equation. In: M. Remoissenet and M. Peyrard, editors, Nonlinear Coher- 
ent Structures in Physics and Biology, Lecture Notes in Physics vol. 393 
( Springer- Verlag, 1991), p. 159-167. 

[18] S. Flach and K. Kladko. Moving discrete breathers ?, Physica D 127, 61-72 
(1999). 

[19] S. Flach and C.R. Willis. Movability of localized excitations in Nonlinear 
Discrete systems: a separatrix problem, Phys. Rev. Lett. 72, 1777 (1994). 

[20] S. Flach, Y. Zolotaryuk and K. Kladko. Moving kinks and pulses: an inverse 
method, Phys. Rev. E 59, 6105-6115 (1999). 

[21] S. Flach and C.R. Willis. Discrete Breathers, Physics Reports 295, 181-264 
(1998). 

[22] D. Hennig, J.F.R. Archilla, J. Agarwal. Nonlinear charge transport mecha- 
nism in periodic and disordered DNA, Physica D 180, 256-272 (2003). 

[23] J.M. Hyman, DW. McLaughlin and A.C. Scott. On Davydov's alpha-helix 
solitons, Physica D 3, 23-44 (1981). 

[24] G. Iooss, M. Adelmeyer. Topics in bifurcation theory and applications, Adv. 
Ser. Nonlinear Dynamics 3, World Sci. (1992). 

[25] G. Iooss, K. Kirchgassner. Travelling waves in a chain of coupled nonlinear 
oscillators, Commun. Math. Phys. 211, 439-464 (2000). 

[26] G. Iooss, E. Lombardi. Polynomial normal forms with exponentially small 
remainder for analytic vector fields. Preprint Institut Non Lineaire de Nice, 
2004. 

[27] G. Iooss, M-C Peroueme. Perturbed homoclinic solutions in reversible 1:1 
resonance vertor fields, J. Diff. Eqs. 102, 62-88 (1993). 

[28] G. James, Y. Sire. Travelling breathers with exponentially small tails in 
a chain of nonlinear oscillators, to appear in Comm. Math. Phys. (2004). 
http: //www.gmm.insa-tlse.fr~james/ 

[29] M. Kastner and J-A. Sepulchre. Effective Hamiltonian for traveling discrete 
breathers in the FPU chain. Submitted to Discrete Cont. Dyn. Sys. B (2003). 

[30] E. Lombardi. Phenomena beyond all orders and bifurcations of reversible ho- 
moclinic connections near higher resonances, Lecture Notes in Mathematics, 
Springer- Verlag, vol. 1741 (2000). 



31 



[31] R.S. MacKay, S. Aubry. Proof of existence of breathers for time- reversible 
or Hamiltonian networks of weakly coupled oscillators, Nonlinearity 7 1623- 
1643 (1994). 

[32] R.S. MacKay, J-A. Sepulchre. Effective Hamiltonian for travelling discrete 
breathers, J. Phys. A 35, 3985-4002 (2002). 

[33] A.M. Morgante, M. Johansson, G. Kopidakis and S. Aubry. Standing wave 
instabilities in a chain of nonlinear coupled oscillators, Physica D 162, 53-94 
(2002). 

[34] M. J.D. Powell. A hybrid method for nonlinear algebraic equations, in Numer- 
ical Methods for Nonlinear Algebraic Equations. Gordon and Breach (1970). 

[35] M. Peyrard, A.R. Bishop. Statistical Mechanics of a nonlinear model for 
DNA denaturation, Phys. Rev. Lett. 62, 2755 (1989). 

[36] M. Remoissenet. Low-amplitude breather and envelope solitons in quasi-one- 
dimensional physical models, Phys. Rev. B 33, number 4, 2386 (1986). 

[37] K.W. Sandusky, J.B. Page and K.E. Schmidt. Stability and motion of in- 
trinsic localized modes in nonlinear periodic lattices, Phys. Rev. B 46, 10 
(1992). 

[38] A.V. Savin, Y. Zolotaryuk, J.C. Eilbeck. Moving kinks and nanopterons in 
the nonlinear Klein-Gordon lattice, Physica D 138, 267-281 (2000). 

[39] A.C. Scott. Nonlinear Science, Oxford University Press, Oxford, 1999. 

[40] J-A. Sepulchre. Energy barriers in coupled oscillators: from discrete kinks 
to discrete breathers, Proceedings of the Conference on Localization and 
Energy Transfer in Nonlinear Systems, June 17-21, 2002, San Lorenzo de 
El Escorial, Madrid, Spain; World Scientihc, Eds L. Vazquez, R.S. MacKay, 
M-P. Zorzano, 102-129 (2003). 

[41] Y. Sire, G. James. Travelling breathers in Klein-Gordon chains, C. R. Acad. 
Sci. Paris, Ser. I 338, 661-666 (2004). 

[42] Y. Sire. Travelling breathers in Klein-Gordon lattices as homoclinic orbits to 
p-tori, submitted to J. Dyn. Diff. Eq. (2004). 

[43] J. Szeftel, G. Huang and V. Konotop. On the existence of moving breathers 
in one-dimensional anharmonic lattices, Physica D 181, 215-221 (2003). 



32 



[44] S. Takeno and K. Hori. A propagating self-localized mode in a one- 
dimensional lattice with quartic anharmonicity, J. Phys. Soc. Japan 59, 3037 
(1990). 

[45] A.J. Sievers and S. Takeno. Intrinsic localized modes in anharmonic crystals, 
Phys. Rev. Lett. 61, 970-973 (1988). 

[46] J.L. Ting and M. Peyrard. Effective breather trapping mechanism for DNA 
transcription, Phys. Rev. E 53, 1011 (1996). 

[47] A. Vanderbauwhede, G. Iooss. Center manifold theory in infinite dimensions. 
Dynamics reported J, new series, 125-163 (1992). 



33 



